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Abstract. We establish a limiting absorption principle for Dirichlet Laplacians in quasi- 
cylindrical domains. Outside a bounded set these domains can be transformed onto a semi-cylinder 
by suitable diffeomorphisms. Dirichlet Laplacians model quantum or acoustically-soft waveguides 
associated with quasi-cylindrical domains. We construct a uniquely solvable problem with perfectly 
matched layers of finite length. We prove that solutions of the latter problem approximate outgoing 
or incoming solutions with an error that exponentially tends to zero as the length of layers tends 
to infinity. Outgoing and incoming solutions are characterized by means of the limiting absorption 
principle. 
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1. Introduction. The perfectly matched layer (PML) method, originally intro- 
duced in [1], is in common use for the numerical analysis of a wide class of problems. 
For some of them stability and convergence of the method have been proved math- 
ematically; see, e.g., [2, 3, 4, 5, 6, 7, 20]. In the present paper we develop the PML 
method for Dirichlet Laplacians in quasi-cylindrical domains Q C K"+^, see Fig. 1.1. 
These are unbounded domains that outside a bounded set can be transformed onto 
a semi-cylinder by suitable diffeomorphisms. Intuitively, one can understand G as 




Fig. 1.1. Examples of quasi- cylindrical domains inK.^. 

a domain whose boundary asymptotically approaches at infinity the boundary of a 
semi-cylinder (0,cx)) x 17, where the cross-section f2 of ^ at infinity is a bounded do- 
main in R". Dirichlet Laplacians A model quantum or acoustically-soft waveguides 
associated with quasi-cylindrical domains. In order to characterize outgoing and in- 
coming solutions of the Helmholtz equation (A — fj,)u = / we establish a limiting 
absorption principle. Then we construct a uniquely solvable problem with PMLs of 
finite length. This is a Dirichlet problem in the domain G truncated at a finite dis- 
tance. We prove that solutions of the latter problem locally approximate outgoing or 
incoming solutions of the Helmholtz equation with an error that exponentially tends 
to zero as the length of PMLs tends to infinity. In other words, we prove stability and 
exponential convergence of the PML method. We find that the rate of exponential 
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convergence depends only on the spectral parameter /x and on the infinitely distant 
cross-section ft. Thus the rate is the same as in the particular case of a domain Q 
that coincides with the semi-cylinder (0, oo) x f2 outside a bounded set. 

As is known, construction of PMLs is closely related to the complex scaling. 
The complex scaling involves complex dilation of variables and has a long tradi- 
tion in mathematical physics and numerical analysis; for a historical account see 
e.g. [16, 10, 34]. Although there are several papers utilizing different approaches to 
the complex scaling in waveguide-type of geometry, e.g. [8. 11, 12, 20], the complex 
scaling has not been used in quasi-cylindrical domains before. Our approach to the 
complex scaling originates from the one developed in [18] for a Schrodinger operator 
in M."^ (see also [16]). Deformations of the Dirichlet Laplacian by means of the com- 
plex scaling give rise to an analytic family of non- self adjoint operators in Q. These 
operators correspond to a Dirichlet problem with infinite PMLs. Localization of the 
essential spectra of these operators togethcir with certain relations between their re- 
solvents justifies a limiting absorption principle. For locating the essential spectra we 
employ methods of the theory of elliptic boundary value problems [24, 25, 28]. Rela- 
tions between the resolvents are obtained with the help of Hardy spaces of analytic 
functions. Note that Hardy spaces in context of the complex scaling were originally 
used in [37, 38]. As we mention in Remark 6.2, our methods also make it possi- 
ble to develop an analog of the celebrated Aguilar-Balslev-Combes-Simon theory of 
resonances. 

As is typically the case, solutions to the Helmholtz equation satisfying the limiting 
absorption principle locally coincide with solutions to the problem with infinite PMLs. 
Moreover, under certain assumptions on the right hand side solutions to the latter 
problem are of some exponential decay at infinity. This allows us to prove unique 
solvability of the problem with finite PMLs and establish exponential convergence of 
the PML method by using the compound expansion technique [27, 25]. In [20] we 
used a similar approach to study the PML method for inhomogeneous media. Here 
we study the PML method for a wide class of quasi-cylindrical domains. 

This paper is organized as follows. Section 2 consists of preliminaries, where 
we introduce notations, formulate our assumptions on the quasi-cylindrical domains, 
and give a formal definition of operators corresponding to the problem with infinite 
PMLs. In Section 3 we demonstrate that the operators are well-defined and derive 
some estimates on their coefficients. As shown in Section 4, these operators give rise 
to an analytic family of m-sectorial operators. In Section 5 we introduce and study 
Hardy spaces of analytic functions. In Section 6 we formulate and prove a limiting 
absorption principle. In Section 7 we show that under certain assumptions on the 
right hand side solutions to the problem with infinite PMLs are of some exponential 
decay at infinity. Finally, in Section 9 we formulate and study the problem with finite 
PMLs and prove exponential convergence of the PML method. 

2. Preliminaries. In this section we introduce basic notations that are in use 

throughout the paper. We formiilatc our assumptions on the quasi-cylindrical do- 
mains and introduce differential operators corresponding to the problem with infinite 
PMLs. Recall that PMLs are artificial strongly absorbing layers designed so that 
waves coming from a non-PML medium to PMLs do not reflect at the interface. 

Let {x,y) and {C,r]) be two systems of the Cartesian coordinates in n > 1, 

such that a;, C S M, while y = (yi, . . . , y„) and rj = (tji, . . . , r/n) are in M". Let dx = 

Consider a bounded domain c M", and the semi-cylinder IR+ x Q, where M+ = 
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{a; e IR : a; > 0}. We say that C C M"+^ is a quasi-cy Under, if there exists a 
diffeomorphism 

R+x09(a;,y)^x(a;,y) = (C,??)eC, (2.1) 

such that the elements >4m(2^' ') Jacobian matrix x' tend to the Kronecker delta 
S^m in the space C°°(fi) as x — > +oo. 

Let Q he a. domain in with smooth boundary dQ. We suppose that the set 

{(C,??) G 5 : C ^ 0} is bounded, and the set {(C,??) G : C > 0} coincides with a 
quasi-cylindcr C. (Extension of our results to the case of a domain Q that coincides 
outside a bounded set with several quasi- cylinders is straightforward.) Following [25], 
we say that ^ is a quasi-cylindrical domain. 

Introduce the notation V^^ = (9^,9^^, . . . ,9^„)^- In the domain Q we consider 
the Dirichlet Laplacian A = — V^^ • V^^, which is initially defined on the set Cq^{Q) 
of all smooth compactly supported functions w in ^ satisfying the Dirichlet boundary 
condition u \dg= 0. 

Consider the complex scaling x ^ x + \s{x — r) with parameters r > and A e C. 
Here s{x) is a smooth scaling function possessing the properties: 

s{x) = for X < 0, (2.2) 
Q<s'{x)<l for all a; e M, (2.3) 
s'{x) = 1 for a; > C > 0, (2.4) 

where s'{x) = dxs{x), and C is arbitrary. For all real A e (—1,1) the function 
M+ 3 X ^ x + \s{x — r) is invertible, and 

M+ X O 9 (a;, y) ^ K,\^r{x, y) = {x + Xs{x - r), y) e M+ x 

is a selfdiffeomorphism of the semi-cylinder. Therefore 

„^ - / ^°'^A,rO^"^(C,??) for (C,??)eC, , . 

^xACv)-^ (C,»?) for (C,^)ea\C, ^^-^^ 

is a selfdiffeomorphism of Q. In other words, -dx^r with A e (—1, 1) and r > is 
a scaling of the quasi-cylinder C along the curvilinear coordinate x. Let (t?^ ,,)^ be 
the transpose of the Jacobian matrix Then e\^r = (''?Ar)^^Ar the matrix 

coordinate representation of a metric eA,r on Q, and 

Aa,^ = - (det eA,r) "'^ • (det eA,r) '^'e^^ V^^ (2.6) 

is the Laplace-Beltrami operator on the Riemannian manifold (^,eA,r)- A-S the pa- 
rameter r increases, the equalities i?A.r(Cj'?) = (Cj'?); ^X-riCiV) = Id, where Id is the 
(n + 1) X (n + 1) identity matrix, and the equality AA,r = A become valid on a larger 
and larger subset of ^. In the case A = the scaling is not applied. Therefore eo,r = e 
is the Euclidean metric and Ao.r = A. 

In order to consider complex values of the scaling parameter A, we impose addi- 
tional assumptions on the diffeomorphism (2.1): 

i. the function 9 a; i— >■ k{x, •) € C°°(0) has an analytic continuation from 
M+ to some sector 



§a = {2 G C : I argz| < a < 7r/4}; 



(2.7) 
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ii. the elements Xg^y^iz, •) of the Jacobian matrix V tend to the Kronecker delta 

Sim in the space C°^{Q) uniformly in z G as z — > oo. 
For instance, the assumptions are satisfied for the following quasi-cylinders: 

C = {(C,7?)eR^(C,J?) = (a;,y + log(x + 2)), x € M+, ye [0,1]}, 

C= {(C,7?) €M^C = _^ ^{t)dt, 7? = #(x), x€R+, 2/G [0,1]}, 

where as (p{x) and ip{x) we can take the functions 1, 1 + e~^, 1 + (a;+ with s > 0, 
1 + 1/ log(x + 2), 1 + 1/ log(l + log(x + 2)), and so on. These examples show that 
quasi-cylinders can have very different shapes comparing with the semi-cylinder. 

In the next section we will show that for all sufficiently large r > the assumptions 
i, ii on X together with (2.2) and (2.3) lead to the analyticity of the coefficients of 
the differential operator (2.6) with respect to the scaling parameter A in the disk 

Pa = {A G C : |A| < sina < l/%/2}. (2.8) 

Thus the equality (2.6) defines A^.r for all A € Va- Clearly, A;^^^ coincides with the 
Dirichlet Laplacian A on the set 

= (g \ C) U {(C, 7?) e C : (C, r?) = Kix, y),x<r,yen}. (2.9) 

We will show that on the set G\Gr the operator A\^r with A G 2?a \ 1^ describes an 
infinite PML for A. In the case SjA > (resp. 9A < 0) this PML is an artificial 

nonreflective strongly absorbing layer for the outgoing (resp. incoming) solutions. 

Remark 2.1. For simplicity we consider in this paper only Dirichlet Laplacians. 
However, similar methods can be used to develop and study PML method for the 
Schrddinger operator A + V inQ, where A is the Dirichlet Laplacian and V G C°°{Q) 
is a real-valued potential with the following properties: for some ro > and a > 
the function x ^ V o k{x,-) G L^(0) extends by analyticity to the sector {z G C : 
I arg(z — ro)| < a}, where for all y (z fl we have \V o >c{z, y)\ < C{\z\) — > as z — > oo. 
One can also include into consideration potentials with moderate local singularities 
and relatively bounded operator-valued potentials. 

3. Construction of infinite PMLs. In this section we show that for all suf- 
ficiently large r > the differential operator (2.6) is well-defined for complex values 
of the scaling parameter A in the disk (2.8). We also obtain some estimates on the 
matrix e^.r- 

Consider the quasi-cylinder C as a manifold endowed with the Euclidean metric 
e. We will use the coordinates {(, rj) in Q and (x, y) in M+ x Q, and identify the 
Riemannian metrics on Q and K+ x fl with their matrix coordinate representations. 
Let g = K*ehe the puUback of the metric e by the diffeomorphism >f in (2.1). Then 
the matrix g = [g£m]"^iLi is given by the equality g = {>c')^ >c', where (V)^ is the 
transpose of the Jacobian k' . Since the diffeomorphism x satisfies the assumptions i,ii 
of Section 2, we conclude that the metric matrix elements 

Sc,3 z^gemiz,-)^C°°m (3.1) 

are analytic functions. Moreover, gem{z,-) tends to the Kronecker delta 6em in the 
space C°°(rJ) uniformly in z G Sq as z — >■ oo or, equivalently, we have 

\d^{gem{z,y) - 6em)\ < Cg{\z\) ^ as z oo, z G §«, t/ G O, \q\> 0; (3.2) 
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Fig. 3.1. The curve il^.r for complex values of X. 



here = d^ld^l • • • ^y" '^i*'^ ^ multiindcx q= {qi,. . . ,qn), and \q\ =J2qj- 

Consider the selfdiffeomorphism Kx^r, A S (—1,1), of the semi-cyhnder R+ x fl. 
We define the metric gA,r = '«A,r*g on M+ x O as the puUback of the metric g by 
K\^r- As a result we get the manifold (1R+ x fl, gA,r) parameterized by A e (—1, 1) and 
r > 0. We deduce 

gx,r{x, y) = diag {1 + As' (a; - r), Id} g{x + Xs{x -r),y) diag {1 + Xs'{x - r), Id} , 

(3.3) 

where Id stands for the n x n-identity matrix, and diag {1 + As'(a; — r), Id} is the 
Jacobian of K^.r- 

Let us consider complex values of the scaling parameter A. We suppose that A is in 

the complex disk , where a is the same as in our assumptions on the partial analytic 
regularity of the diffeomorphism >^; cf. (2.7), (2.8). The curve £A,r = {z & C : z = 
X + \s{x — r),x > 0} lies in the sector S^, see Fig. 3.1. We define the matrix g^.r for 
all non-real A in the disk by the equality (3.3), where g{x + Xs{x — r) , y) stands for the 
value of the analytic in z G function g{z,y) a,t z = x + Xs{x — r). By analyticity 
in A we conclude that g^,,- is a complex symmetric matrix, the Schwarz reflection 
principle gives g^.r — gjr^ where the overlinc stands for the complex conjugation. If 
A € Va is non-real, then the matrix gA,r does not correspond to a Riemannian metric. 
However, gA,r is invertible for all A e Va provided r > is sufficiently large. Indeed, 
s{x — r) = and g^^.(a;, y) = g~"'^(x, y) for all x < r. On the other hand, the matrix 
g(a; + Xs{x — r),y) from (3.3) is invertible for a; > r with large r > as it is only little 
different from the identity matrix; the latter fact is a consequence of (3.2) and the 
inequality \x + Xs{x — r)\ > r. 

The derivatives d^d^g^l. are analytic functions of A S Va. From the condi- 
tions (3.1) and (3.2) together with (2.4) and (3.3) it follows that 

^x5^(gA.r(-7^-?/)-diag{(l + A)-2,Id})|| ^Oasx^+oo, y G H, A G (3.4) 

and the estimate 

5^5«(g-i(x,j/) -diag{(l + As'(x-r))-2,Id})|| < Cp,(r) (3.5) 

holds uniformly in {x,y) G [r,oo) x and A € Pa, where || • || is the matrix norm 
ll^ll = max^TO \(iem\, and p + |g| > 0. The constants Cpq(r) in (3.5) tend to zero as 
r — +00. 

We define the complex scaling dx^r for all A in the disk Va by the equality (2.5), 
where HOKx,r{x, y) is the value of the analytic in 2 G function k{z, y) at the point 
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z = X + Xs{x — r). Consider the matrix e^.r = r)^^A r- clear that for all 

(CtV) G G\C the matrix eA,r(C)^) coincides with the (n + 1) x (n + l)-identity. For 
all real A e we have 

g\,r{x,y) = {>i'{x,y))''' {ex^r° >i{x,y))>€'{x,y), {x,y)GR+xfl, (3.6) 

where gA,r = x*^\,r is the puUback of the corresponding metric e^.r on C by the dif- 
feomorphism x. Therefore e^.r is analytic in A G I?q and invertible for all sufficiently 
large r > 0. By analyticity in A we conclude that e\^r{C,v) is a complex symmetric 
matrix, the Schwarz reflection principle gives = ^ r- 

Differentiating the equality (3.6), we see that the matrices d'^d^ex,r and d^d^e^^^. 
are analytic in A e V^. Moreover, from (3.4), (3.5), and our assumptions on x we 
obtain 

^C^^KUC^y) - diag{(l + A)-', Id}) II ^ as C ^ +00, 
d^d^vKriC, V) - diag{ (1 + As;(C, r,)) Id}) < c(r). 



where p + |g| < 1, and c(r) — ^ as r — )• +oo. 



(3.7) 



Here (C, 77) e C and s'^{C, rj) stands for the function s'{x — r) written in the coordinates 
(C, ??)■ Wc extend from C to C/ by zero. Note that the estimate (3.7) remains valid 
for all (C, Tj) ^ G and the constant c(r) is independent of A € and (C, rj) G Q. Now 
we see that for all sufficiently large r > the differential operator (2.6) is well-defined 
for all A in the disk Va, and its coefficients are subjected to the estimate (3.7). 

4. Analytic families of operators. In this section we study the unbounded 
operator A^.r in the Hilbert space L'^{Q) with the usual norm 

\\u;L\g)\\ = (^jHC,ri)\''dCdrij'\ 

The operator A^.r is initially defined on the set C^{Q). We show that the operator 
is closable, and the closure defines an analytic family 9 A A^.r of type (B) [23]. 
In a standard way this implies that the resolvent (A^.r — A*)~^ is an analytic function 
of A and ^ on some open subset of Va x C. The latter fact will be used for justiflcation 
of a limiting absorption principle in Section 6. 

We intend to show that the operator A^^r in L'^{G) with the domain Ct^{Q) is 
sectorial and to study its Friedrichs extension. With the operator A^.r we associate 
the quadratic form 

q'A,r[w,M]= / ((deteA,r)^^^e^^^Vc,,u, Vc^(dete5^ (4.1) 
Jg 

where (•,•) is the Hermitian inner product in C"+^ and u G C^{Q). Let (•,•) stand 
for the inner product in L'^{Q)- We represent the quadratic form as follows: 

qX,r[u,u] = {-^CV ■ SA,r^C'7"' 

f 12 12 ("^ ^) 

+ / {{detex,r)^^^exl'^ir,u,uyQr,{dete^^y^^^)dCdr]. 
Jg 
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For the first term in the right hand side of (4.2) we prove the following lemma. 

Lemma 4.1. Assume that r > is sufficiently large. Then there exist ip < 7r/2 
and 6 > such that for all A G Va and u G C^{Q) we have 

|arg(-Vf,, • e^^^Vf^w,u)| < <p, 5{Au,u) < iR(-Vc^ • e^_^V^^u,u) < 5~'^{Au,u). 

In other words, the form (— V^,, • e^^V^,,u, m) ^5 sectorial. 
Proof. It is clear that 



Jc 



(4.3) 



Let us estimate the numerical range of the matrix e^^(C,77). We shall rely on the 
estimate (3.7). Let ^ = Vc,,u(C,?7) e C"+^ Observe that by virtue of < s'^iCrj) < 1 
and |A| < sina < 2~^/^ we have 



1^1 V4< ^diag{(l + As;(C,77)) ',IdK < 121^1' 



arg(e-diag{(l + As;(C,r,)) ',Id}^ 



< 2a. 



(4.4) 



Since r is sufficiently large, the constant c(r) in (3.7) is small. In particular c(r) 
meets the estimate 4(n + l)^c(r) < sin(cr/2) with some a G (0, 7r/2 — 2a). Then (3.7) 
together with (4.4) gives 



arg(^-e-i(C,r;)e) <<p<7r/2, 6]^ < ^ ■ e^^^iCrj)^ < S-'\^\^ (4.5) 
uniformly in A G Va and {(, rf) G Q, where = 2a + u and 

5 = min|l/4 - (n + lfc{r), (l2 + (n + lfc{r)) |. 

Taking into account (4.3) we complete the proof. □ 

Remark 4.2. Throughout the paper we say that r > is sufficiently large if the 
matrix e\_r{C,rj) '-^ invertihle and its inverse meets the estimates (4.5) uniformly in 
iCv) eg and A G ©a. 

In the next lemma we show that in the right hand side of (4.2) the second term 
has an arbitrarily small relative bound with respect to the first term uniformly in 

A e Pa- _ 

Lemma 4.3. For any e > and u € Cq°(5) the estimate 



/ ((deteA,r-)^^\r^Cr,w,wV^r,(dete^ J ^^^)dCdr] 
Jg 



< e|(-Vc, • Vc,«,w)| + Ce-'\\u- L\g)\[' 



holds, where the constant C is independent of e, u, and A G I?a • 
Proof We have 

/ ((deteA,r)^^^e;^^Vc^u,MVc,,(dete^^) ^^^)dCdr] 

Jg 

- ^Ug K^A,rVcr,^^,Vf^(dete^,)-'/')|2dCrfr?)'^'(^ \u\^dCdrj) 
< c{ej |(e^;,Vc^«,Vc,(dete^_J-'/')|2(iCdr, + e-i||w;L2(g)||2 



1/2 



(4.6) 
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with arbitrarily small e > and C = sup(^ (det eA,r(C) 'y)) < oo, cf. (3.7). 
Prom (4.5) it follows that 

|^{?-e^_i.(C,?7)0| < (tan<^)SR{e-e^,UC,'7)C}, 

where the form 3?{^ • ^x].{C,il)£.} defines an inner product in C"+^. This and the 
Cauchy-Schwarz inequality give 

\r-^x!ri('V)i\^ < (l+tan^)2K{e-e^i(C,r;)e}3?{r-e^i(C,r?)r}, (4.7) 

cf. [23, Chapter VI.2]. Wc substitute ^ = Vc,,m(C,??) and r = Vc,,(deteA,r(C,»?))~^^^ 
Thanks to (3.7) we have the uniform bound 

|Vc„(deteA,r(C,»?))"'^'r <c, A e P^, {C,v)&G- 
This together with (4.7) and (4.5) implies 

K^A.^^C')"' Vcr,(dcte^ ,,)"^/^)|2 < c5-i(l +tan(^)23fj(e^^i^Vc^u, Vc„u). 

Now we make use of (4.6) and establish the assertion for C = C^C(5~^(1 +tan(p)^ and 
an arbitrarily small e = Cc5~^(l + tanip)^e. □ 

As a consequence of the equality (4.2) and Lemmas 4.1, 4.3 for all sufficiently 
large r > we obtain 

\a.vg{qxAu,u] + j\\u;L\g)f)\<^<7T/2, (4.8) 

SqoAu,u]-j\\u;L^{g)f < 3?gA,r[w,w] < 6-\qoAu,u] + \\u;L^{g)f) (4.9) 

with some angle (p and some positive constants 6 and 7, which are independent of 
u G C^{Q) and A G Va- The symmetric form qoAu,u] = Jg(^ Qr)U,V Qr^u) dC,dr} 
is independent of r, it corresponds to the Dirichlet Laplacian A = Aq^^- Clearly, 
qxA^-iU] = (AA,rW,w). Estimate (4.8) implies that the numerical range 

{/X e C : M = {Ax,rU,u),u& C^iQ), \\u;L^{g)\\ = 1} 

is a subset of the sector {/x G C : | arg(^ + 7)| < < 7''/2}. By definition this means 
that the operator AA.r with the domain C^(Q) is sectorial. 

We introduce the Sobolev space H^{Q) as the completion of the set C^{Q) with 
respect to the norm 

\\u;H\g)\\ = ^qoAu,u] + \\u;L^{QW- 

Recall that i) a sequence {uj} is said to be (/a, r— convergent, if Uj is in the domain 
of qx,r, 11% — zt;L^(5)|| — >■ and q\,r[^j ~ Um,Uj — Um] — >■ as j, m — >■ 00; n) the 
form qx,r is closed, if every qrA.r-convergent sequence {uj} has a limit u in the domain 
of qx,r, and gA,r[u ~ — Uj] 0. From (4.8), (4.9) it immediately follows that 

qx.r with the domain H^{Q) is a closed densely defined sectorial form [23, 34], and 
its sector {/x G C : | arg(/x + 7)! < (p} is independent of A G 2?a- As known [23, 
Chapter VI. 2.1], to every closed densely defined sectorial form there corresponds a 
unique m-sectorial operator. Namely, to the form qx,r there corresponds a unique 
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m-sectorial operator A;^ in L^{A4) such that its sector is the sector of qx^r^ the 
domain D(AA,r) is dense in H^{Q), and qA,r['W,u] = {A\^rU,v) for all u e D(AA,r) 

o 

and V G H (Q). (Here and elsewhere m-sectorial means that the numerical range 
{fj, = {Au,u)-n £ C : u £ D(A),(u,u)^ = 1} and the spectrum a{A) of a closed 
unbounded operator ^ in a Hilbert space H with the inner product (•, ■)-h both lie in 
some sector {^u e C : arg(/x — j) < if} with 7 e R and (p < n/2.) In particular, to the 
symmetric nonnegative form go,r there corresponds the selfadjoint Dirichlet Laplacian 
A = Ao^r- The m-sectorial operator A>^ ^ in L'^{Ai) is the Priedrichs extension of the 
sectorial operator A^ ,. defined on C(^{Q), see [23, Chapter VI. 2. 3]. As we show in 
assertion 1 of the next proposition the set Cq^{G) is a core of the Friedrichs extension. 

Proposition 4.4. Assume that r > is sufficiently large. Then the following 
assertions are valid. 

1. For A e Va the m-sectorial operator AA,r in Li^iQ) with the domain D(AA,r) 
is the closure of the operator A\^r defined on the set C^{Q). 

2. The family of m-sectorial operators Da 9 A A^.r in L'^{Q) is an analytic 
family of type (B) . 

3. The resolvent T B (A,/x) i-)- (A^^r — A*)~^ : L'^iG) ~^ L'^iG) o-n analytic 
function of two variables on the set F = {(A, /i) : A e ©a,/" € C \ a{A\^r)}, 
where a{Ax,r) is the spectrum of Ax,r- 

Proof. Consider the domain D{A\^r) as a Hilbert space endowed with the graph 
norm Hu; D(AA,r)|| = \\u-,L'^{G)\\ + || AA,rM; -^^(^)|| ■ Let n he a point outside of the 
sector of the m-sectorial operator A^.r- Then the set 

C(Aa,.) = {u:u = {Ax,r - lA~^fJ& 

is dense in D(AA_r) because the resolvent (A^^r — • L'^iG) ~^ D(AA,r) is bounded 
and the set Cq'{G) is dense in L^{G). From (4.5) it follows the estimate 

5R(C-eA,'(C,'?)0 >ckl^ CeM"+SAePa,(C,?7)ea, (4.io) 

on the principal symbol of Aa.d where c > 0. Hence AA,r is a strongly elliptic 

operator. As is well-known, a strongly elliptic operator and the Dirichlet boundary 
condition set up an elliptic boundary value problem, e.g. [26]. The usual argument 
on the regularity of solutions to the elliptic boundary value problems [26, 25] implies 
that the set C(AA,r) consists of smooth in Q functions u with u \dG= 0. Multiplying 
u <E C(AA,r) by appropriate cutoff functions Xj with expanding compact supports 
suppxi C suppxi-i-i, it is easy to see that for any u e C(AA,r) there is a sequence 
{Xj'^}j^i such that XjU S Cg^{Q) tends to u in D(AA,r) as j -\-oo. Assertion 1 is 
proven. 

The family 3 X q\^r is analytic in the sense of Kato (i.e. q\^r is a closed 

o 

densely defined sectorial form, its domain (^) is independent of A, and the function 

o 

T>a 3 A qx,r[u,u] is analytic for any u G H {Q)). By definition [23, 34] this means 
that the family of m-seetorial operators Pq, 9 A i— >■ Aa^t is an analytic family of type 
(B). This proves assertion 2. 

As is well-known [23, 34], any analytic family of type (B) is also an analytic family 
of operators in the sense of Kato. Now a standard argument justifies assertion 3; see, 
e.g., [34, Theorem XH.T]. □ 

5. Spaces of analytic functions. We have shown that for all fj, outside the 
sector of the family 9 A 1-^ AA,r of m-sectorial operators the resolvent (AA,r — 
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is an analytic function of A, see Proposition 4.4.3. In order to get some relations 
between the resolvents (A — ^)~^ and (A;^ — ^)~^ we will use a sufficiently large 
Hilbert space Ha{G) of analytic functions 

foi3x,r&L\g); (5.1) 

here i^^.r is the complex scaling (2.5). The goal of this section is two-fold: 

1. To introduce a Hilbert space Ha{G), which is sufficiently large in the sense 
that for any A G X>„ and r > the set {/ o ■dx,r e L'^{g) : f G "Hal^)} is 
dense in L'^{Q)\ 

2. To derive the uniform in A € and / e 'Ha{Q) estimate 

\\fo'9xy,L\g)\\<Cr\\f;Ua{G)l r>0. (5.2) 

Introduce the Hardy class i5(§a) of all analytic functions 9 ^ F{z) G L'^{Sl) 
satisfying the uniform in tp estimate 



/>oo 

/ \\F{e''^x);L^{n)fdx<CF, ^G(-a,a). (5.3) 
Jo 



Below wc cite some facts from the theory of Hardy classes. 
Lemma 5.1. 

1. Every F e Sj{Sa) has boundary limits F± e L^(M+ x O) such that 



/ \\F{e''''x) - F±{xy, L'^{n)f dx ^ as ip ^ ±a. 
Jo 



2. The Hardy class ^{Sa) endowed with the norm 

\\F;Sj{S^)\\ = {\\F_-L'{R+ x n)f + \\F+-L'{R+ x n)ff^ 
is a Hilbert space. 

3. For every compact set AcSa there is an independent of F & Sj{Sa) constant 
C{R) such that for all z G R we have 

\\F{zy,L\n)\\<cimF;Sj{^a)\\. 

4- For any F e i5(§a); z G Sa, and ip G (—a, a) we have 

\\Fiz + e'^x); L\n)fdx < \\F;S)iS^)f. 



Proof Assertions 1-3 are direct consequences of well-known facts from the theory 
of Hardy spaces of functions analytic in strips, e.g. [35]. In fact, the proof reduces to 
the conformal mapping of the sector §„ to the strip {z G C : —a < Qz < a} , we omit 
the details. 

Let us prove assertion 4. A standard argument, see e.g. [36], shows that any 
F G Sj{Sa) can be represented by the Cauchy integral 

Fiz)=± l^^^^dx-1-. /^T!^,,, ,,s„, (5.4) 
^ 2m J z- e»«x 2Tn J z- e-'"x a, v / 

— oo — oo 
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where we assume that the boundary Umits F±{x) are extended to x < by zero. 
The integrals arc absohitcly convergent in the space L'^{fl). As is well-known [35], 
the first integral in (5.4) defines an element /+ of the Hardy space io(C^) of L'^{Q)- 
valued functions analytic in the half-plane C~ = {z € C : 9(e~*"z) < 0}. As shown 
in [37, 38], the norm in Sj{C~) can be defined by the equality 

sup r \\U{z+ + e'^xy,L\n)fdx] , e--^+eM. 

i/>G(a-7r,a) Jo / 

Then ||/+;^(C-)|| = \\F+; L^{R+ x Similarly, the second integral in (5.4) defines 
a function /_ from the Hardy space i5(C+) in C+ = {z £ C : S(e'"2;) > 0}, and 
||/_;i5(C+)|| = \\F_;L^R+ x f2)||, where 



||/_;i5(C+)||= ( sup r\\Mz.+e'^xy,L\n)fdx 



1/2 







Q,7r— a) JO 

As a consequence, for all z &Sa and ip e {—a, a) we have 

||F(z + p/^x);i2(f7)||2dx<||/+;io(Q)f +||/_;i^(Cj)f 
= \\F+;L\R+ X Q)f + \\F_;L^{R+ x 17)f = ||F;.^(§„)f . 

□ 

Consider the algebra ^ of all entire functions C B z i-^ F{z) € C^{^) with the 
following property: in any sector \'^z\ < (1 — e)3fiz with e > the value ||i^(z); L?{^)\\ 
decays faster than any inverse power of as "^z — >■ +00. Examples of functions 
F £ S are F{z) = e"'*'^ P{z), where 7 > and P{z) is an arbitrary polynomial in z 
with coefficients in C^(J7). Clearly, S C i5(Sc«)- The next lemma is an adaptation 
of [18, Theorem 3], we omit the proof. 

Lemma 5.2. The set of functions {R+ x Vt ^ {x,y) ^ F o K\^r{x,y) : F e S"} is 
dense in the space L^(K+ x il) for any A € Va and r > 0. Here F o K\^r(x, ■) stands 
for the value of the entire function z h- > F(z) at the point z = x + Xs{x — r). 

Now we are in position to prove the following proposition. 

Proposition 5.3. 



1. The estimate 



f 

Jo 



\FoKxy,L^in)fdx<Cr\\F;^{Sa)f, r > 0, (5.5) 



holds uniformly in X € T>a and F € Sj{Sa). 

2. The space SjiSa) is sufficiently large in the sense that for any A S T>a and 
r > the set {R+ x fl B {x,y) F o K\^r{x,y) ■ F e i5(Sa)} is dense in 

l2(m+ X n). 

3. For any F G S){Sa) and r > the function S A H> F o AtA.r S L^(R+ x fi) 
is analytic. 

The proof is preceded by a discussion. By Proposition 5.3.1 for any r > the 
complex scaling K,\_r induces the uniformly bounded injection 



S){Sa)B F^ FoKx,r&L^{^+^^), XeT>c. 



(5.6) 
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Any function F G Sj{§a) can be reconstructed from its trace x FoK\^r{x, •) € L^{fl) 
by analytic continuation from the curve {x + X.s{x — r) G C : a; > 0} to the sector Sq,. 
Therefore we can always identify the space Sj{Sa) with the range of injection (5.6). 
By Proposition 5.3.2 the range is dense in L^(IR+ x fi). 

Proof. 1. For all A G Va the curve {x + Xs(x — r) E C : x > 0} lies in the sector 
Sa- Observe that this curve is differ from a ray only inside an independent of A G Va 
compact subset Mr C §«• Now the uniform in A G Va and F G i5(Sa) estimate (5.5) 
follows from assertions 3 and 4 of Lemma 5.1. 

2. The assertion is an immediate consequence of the embedding # C Sj{E>a), 
Lemma 5.2, and the estimate (5.5). 

3. It is easy to sec that the function Dq, B X i-^ F o Kx,r € Z/^(M_|_ x CI) is weakly 
(and therefore strongly) analytic. □ 

Introduce the Hilbert space 'Ha{Q) with the norm 

\\f;nam = \\f:L\g)\\ + \\f o K--^-^{^a)\\. 

The space T-LaiQ) consists of all functions / G L'^[Q) such that / o is an element 
of the Hardy space S){^a)', here x is the diffeomorphism (2.1). As a consequence of 
Proposition 5.3 and definition (2.5) of the complex scaling ■!?A,r we immediately get 
the following assertions. 
Corollary 5.4. 

1. For any r > the estimate (5.2) holds with an independent of f & 'Ha{G) 
and A G constant Cr- 

2. For any A G ©a and r > the set {/ o 'd\^r € L'^{Q) ■ f G J^aiG)} is dense 
in the space L^{Q). 

3. For any f G HaiQ) and r > the function (5.1) is analytic. 

6. Limiting absorption principle. Introduce the Sobolev space Hq{G) of func- 
tions satisfying the homogeneous Dirichlet boundary condition on dQ as the comple- 
tion of the core Cq^{Q) with respect to the graph norm 

\\u;H^{g)\\ = \\u;L'{g)\\ + \\Au;L\g)\\. (6.1) 

(Integrating by parts and using the Cauchy-Schwarz inequality one can easily see that 
the norm (6.1) is equivalent to the traditional norm (X]£+|m|<2 W^^dl^u; L'^{g)\\^y^'^ .) 
By Proposition 4.4.1 the space Hq{Q) is the domain D(A) of the selfadjoint Dirichlet 
Laplacian A. For the points fj, G cr(A) the resolvent (A — /x — ie)~^ does not have 
limits in the space of bounded operators .!^{L?{Q),Hq{Q)) as e tends to zero from 
below (e t 0) or from above (e 4 0). However the limits may exist in the space of 
bounded operators acting from a smaller source space to a larger target space. As 
a source space we take the Hilbert space HaiG) constructed in the previous section. 
As a target space we take the reflexive Frcchet space i?oioc(^)- The space i?oioc(^) 
consists of all distributions u such that gu G Hq{Q) with any g G C^{g), the topology 
in loc(^) induced by the family of seminorms u i-)- \\gu; HQ{g)\\; here C^{g) is 
the set of all smooth functions with compact supports in g. 

Theorem 6.1. Let cr(Ao) stand for the spectrum of the selfadjoint Dirichlet 
Laplacian Aq in L'^{fl). Assume that /^o S M \ a{An) is not an eigenvalue of the 
selfadjoint Dirichlet Laplacian A in L'^{g). Then the following assertions hold. 
1. For all sufficiently large r > and X G Va \ R the resolvent 



{Ax,r - i^o)-^ L\g) ^ Hl{g) 
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is a bounded operator. 

2. The resolvent (A — //q — ie)~^, e ^ 0, viewed as a bounded operator acting 
from Ha{Q) to -f^o ioc(^)' limits as e 10 and e f 0. 

3. Suppose that r > is sufficiently large, X G Va \ K, and f S V-aiG)- Let 
u\^r G Hq{Q) he given by the equality w^.r = {^\,r^ lk))^^{.f°''3x.r)- Then the 
outgoing U- € Hq^^^{Q) and the incoming u+ G -f^oioc(^) solutions defined 
by the limiting absorption principle 

= lim(A — fJ,o — i^)~^f, U- = lim(A — fio — ie)^^f, (6.2) 



meet the relation 



Here the bounded domain Qr is the same as in (2.9). 

The proof is preceded by a discussion. From Theorem 6.1.3 we see that the 
equation {Ax^r — IJ-ajux^r ~ .fo^x,r with non-real parameter A G Va describes infinite 
PMLs on Q \Qr for all /io satisfying the assumptions of the theorem. The layers are 
perfectly matdicid in the sense that for SA > (resp. for 3A < 0) MA,r coincides in Gr 
with the outgoing solution U- (resp. with the incoming solution 7i+). The PMLs are 
absorbing because in contrast to u± the function ux^r decays at infinity in the mean 
as an element of Hq{Q). In the next section we will refine results of Theorem 6.1 
by showing that under an additional assumption on / o the solution ux,r is of 
some exponential decay at infinity. For instance, this assumption is a priori met for 
/ € LP'iQ) supported in Qy Then / o 'dx^r = / and the operator A^.r completely 
describes infinite PMLs on ^ \ • 

Proof. The proof consists of two steps. 

Step 1. In Section 8 below we will show that the graph norm of A;^^^ is an 
equivalent norm in Hq{Q). This immediately implies that D(AA^r) = Hq{Q) as the set 
C^{Q) is dense in both spaces, cf. Proposition 4.4.1. We will also localize the essential 
spectrum (Jess{A.x,r) of the unbounded m-sectorial operator A^.r in L'^{Q). As is well- 
known, the spectrum fT(Ao) consists of infinitely many positive isolated eigenvalues. 
It turns out that aess{A.x.r) consists of an infinite mmibcr of rays emanating from 
every point u G a(Ao), cf. Figure 6.1. By definition c7(An) is the set of thresholds of 




• thresholds v G a-(An); 

rays of cress(AA,r)- 



Fig. 6.1. Essential spectrum of the m-sectorial operator A^,,- for OA > 0. 

the Dirichlet Laplacian A. As A varies, the ray {/x G C : arg(/x— i/) = — 2arg(l + A)} of 
the essential spectrum cress (A^.r) rotates about the threshold u e f7(Ao) and sweeps 
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the sector {/x e C : | arg(^ — < 2a}. In order to avoid several repetitions we 
organized the paper so that the proofs of these results in some greater generality are 
postponed to Section 8; here we take these results for granted. 

Recall that /x is said to be a point of the essential spectrum aess (A) of a closed 
unbounded operator A in the space L'^(G) with the domain Hq{Q), if the boimded 
operator A — n : Hq{Q) — >• L'^{G) is not Fredholm (a linear bounded operator between 
two Banach spaces is Fredholm, if its kernel and cokernel are finite dimensional, and 
its range is closed). Since the operator Aa.t is m-sectorial, there exists a regular point 
of Ax^r in the simply connected set C \ aess{A.x^r)- Therefore 

C \ c7ess{Ax,r) 9 M ^ A^,^ - : H^{g) ^ L^Q) 

is a Fredholm holomorphic operator function, e.g. [24, Appendix]. Recall that the 
spectrum of a Fredholm holomorphic operator function consists of isolated eigenvalues 
of finite algebraic multiplicity, e.g. [24, Proposition A. 8.4]. As a consequence, the 
resolvent 

C \ ae..(AA,.) 9 M ^ (AA,r - m)"' : L^g) ^ H^{g) (6.3) 

is a meromorphic operator function. 

For f,gG L'^iG) and a sufficiently large r > we define the quadratic form 

{f,9)x,r= / /^A/det ex,rdCdr], A € 
Jg 

This form is bounded in L'^{Q). Indeed, thanks to the estimate (3.7) on e^^ we have 

< Ci < I deteA,r(C7'?)l ^ ^2 uniformly in A € Va and (C,^/) "= Q- 

Assume that A e Va is real. Then the form (•, ■)x.r is the inner product induced 
on Q by the metric Qx,r, the norm ^ (/, f)x,r is equivalent to the norm ||/;L^(C/)||, 
and A^.r is the Laplace-Beltrami operator on {Q,ex,r)- We have 

{A-^^)u={{Ax,r-^i){uo^x,r))o'&ll^ \lu&C^{g). (6.4) 

Assume that jj, is not in the sector of A^.r- Then (A^.r — A*)~^ is a bounded operator 
and we can rewrite (6.4) in the form 

(A - = ((Aa,. - lA-\f o ^A,r)) o ^l^^, (6.5) 

where / is in the set {/ = (A — /i)u : u G C^{Qy\. This set is dense in L'^{Q), 
because C^{Q) is dense in Hq{Q), and the operator /S. — ji : Hq{0) — > L'^{G) yields 
an isomorphism. It is clear that (/ o '&x,r, f ° '&\,r)\,r = (/, /)• As a consequence, the 
(real) scaling / i— > / o 'dx,r realizes an isomorphism in LF'iQ), and the equality (6.5) 
extends by continuity to all / G L'^{G)- Taking the inner product of the equality (6.5) 
with g e L'^iG), and passing to the variables (C, »?) = ''^xiC^v) in the right hand side, 
we obtain 

((A-M)-V,5) = {{Ax,r-f^)-\fo^x,r),9 0^x,r)x,r- (6-6) 

Now we assume that /, g e TLaiG)- Then f odx,r and godx,r are i^(t/)-valued analytic 
functions of A in the disk Va, see Corollary 5.4.3. This together with Proposition 4.4.4 
implies that the right hand side of (6.6) extends by analyticity from A e Va fl R to all 
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A e T^a- The right hand side of (6.6) extends from all // outside of the sector of A^^^ 
to a meromorphic function of /i G C \ (Tesj;(A;!^ ,.). In particular, for all A G \ IK. we 
have Ho € <C\aess{^x,r)i cf- Figure 6.1. Here /xq is the same as in the formulation of 
the theorem. 

Now we are in position to prove assertion 1. Consider the projection 
P = s-limie(A — hq — ie)~^ 

e-lO 

onto the eigenspace of the selfadjoint operator A. Suppose, by contradiction, that the 
resolvent (6.3) has a pole at the point /io G C \ <Jess{^x,r)- By Corollary 5.4.2 there 
exist / and g in the space HaiQ) such that /io is a pole of the right hand side of (6.6). 
The equality (6.6) implies that (Pf,g) 0, and thus ker(A — /io) {0}- This is a 
contradiction. Assertion 1 is proven. 

Step 2. We need to show that for any g G C'^{Q) the operator q{A — — ie)~^ 
tends to some limits in the space of bounded operators ,^{J^/fa{G) , Hq (Q)) as e J, 
and e t 0. We take a sufficiently large r = r{g) > such that supp C Gr- Then 
8 o i?A,r = Q for all A e Va- Now we can pass from (6.5) to the equality 

g{A - A*)-!/ = g{A^,. - fi)-\f o i9,,,). (6.7) 

For / e 'Ha{Q) the equality (6.7) extends by analyticity to all A € Consider, for 
instance, the case SJA > (the case SJA < is similar). Then the upper half-plane 
C"*" = {/i G C : 3/i > 0} and a complex neighborhood of the point /io do not contain 
points of cress(AA,r), cf. Figure 6.1. Therefore the right hand side of (6.7) has a 
meromorphic continuation in /t to the union of C"*" and a complex neighborhood of 
/to. Hence the left hand side of (6.7) has the same meromorphic continuation. Clearly, 
a pole at /to may only appear due to a pole of the resolvent (6.3) at /to, but it is a 
regular point by assertion 1. Since q is an arbitrary smooth function supported in Q^- 
this proves assertion 3. In order to prove assertion 2 it remains to note that 

\\\img{A-ii^-ie)-^f-Hl{g)\\ = \\e{A^^, - ii^)-^{f o{}^^,)-Hl{g)\\ 
64,0 

<C{g)\\{A^,r-^io)-']L\g)^Hl{g)\\\\fodxy,L\g)\\ 

<C{g,r,iio)\\f;Ua{g)\\. 

In the last inequality we used Corollary 5.4.1 and assertion 1. □ 

In the following two remarks wc collect some results that can be obtained by 

methods developed in the proof of Theorem 6.1. Although these results are not used 

in this paper, they provide additional insights of the problem. 

Remark 6.2. On the basis of the equality (6.6), Corollary 5.4, and the description 

of aess{A\ r) for \ G 'Da, one can develop an analog of the celebrated Aguilar-Balslev- 

Combes- Simon theory of resonances [10, 16, 18, 34]- We announce some results below, 

for the proof we refer to [21 ]. 

1. The selfadjoint Dirichlet Laplacian A in g has no singular continuous spec- 
trum, its eigenvalues can accumulate only at thresholds v G (T(An). (Exam- 
ples of accumulating eigenvalues can be found e.g. in [13].) 

2. The spectrum (j(AA,r) of the m-sectorial operator A\ r does not depend on the 
choice of the scaling function satisfying (2.2) (2.4). Moreover, the spectrum 
a{Ax^r) liss in the half-plain C+ in the case 9A > and a{Ax^r) C C~ in 
the case 9A < 0, where = {/x G C : 3/t ^ 0} and A G Va- 



16 



VICTOR KALVIN 



3. A point G R \ cr(AQ) is an eigenvalue of A if and only if it is an isolated 

eigenvalue of A\ r, where A G \ K a,nd, r > is sufficiently large. 
4- The resolvent matrix elements ((A — iJ,)~^f,g), where f,g & ^HaiQ), have 
meromorphic continuations from the physical sheet <C\(Jess (A) across fTess(A) 
to the set C \ cfess{A\.r)- Moreover, ^ is a pole of the continuation for some 
f,g € HaiG) if and only if it is an isolated eigenvalue of Ax^r- The non-real 
isolated eigenvalues of Ax^r are naturally identified with resonances of A. 
5. As A changes continuously in the disk T>a, an isolated eigenvalue of Ax^r 
survive while it is not covered by one of the rotating rays of aess{^\,r)- 
Remark 6.3. The argument of the second step in the proof of Theorem 6.1 allows 
also to see that the resolvent (A - ij)-^ : Hl^^^{g) h as a meromorphic 

continuation from the physical sheet C\(Tess(A) across the intervals {u_,u_^_) between 
the neighboring thresholds u± e a"(Aj2) to a Riemann surface. The surface consists 
of the physical sheet C \ (Tess(A) of the Dirichlet Laplacian and an infinite number 
of the sectors {/x € C : > arg(/z — V-) > —2a} attached to C"*" C C \ CTess(A) and 
of the sectors {n G C : < arg(/i — i/-) < 2a} attached to C~ C C \ aess{A) along 
the intervals {i>-^,i>+). Indeed, for any g € C^{g) we can take a sufficiently large 
r = r{Q) > such that suppg C Qr- As A varies in T>a n C"*" (resp. in T>a n C~) 
the strip between the neighboring rays {/z e C : a.rg{z — i'±) = — 2arg(l + A)} of 
(Tess(AA.r) swccps the scctor {/i € C : > arg(/i — > —2a\ (resp. the sector 
{/U G C : < arg(/i — < 2a}) and the right hand side of (6.7) provides the left 
hand side with a meromorphic continuation to the strip. The poles of the continuation 
are resonances of the Dirichlet Laplacian A [39]. 

The results listed in Remarks 6.2 and 6.3 are new and might be of their own 
interest. Traditionally, when studying Laplacians in the waveguide- type of geometry, 
one imposes more or less restrictive assumptions on the rate of convergence of the 
metric g on (0, cxd) x Q to its limit at infinity; see, e.g., [9, 11, 12, 13, 14, 15, 19, 29, 30, 
31]. Contrastingly, our assumptions on the diffeomorphism >c allow for arbitrarily slow 
convergence of the metric g }<*e to the Euclidean metric e at infinity, sec Section 3. 
As a substitution for the assumptions on the rate of convergence of the metric g at 
infinity we use assumptions on the analytic regularity of the diff'eomorphism >c. 

Let us also mention here the paper [33], where general elliptic problems whose 
coefficients slowly converge to their limits at infinity are considered. It is shown that 
any finite accumulation point of eigenvalues corresponding to exponentially decaying 
eigenfunctions is a threshold, these accumulations are characterized in terms of some 
non-classical "augmented scattering matrices." An additional investigation on decay 
of eigenfunctions at infinity is required in order to say whether these results describe 
accumulations of eigenvalues of the Dirichlet Laplacian A or not. This goes beyond 
the scope of the present paper, we refer to [22]. 

7. Exponential decay of solutions in infinite PMLs. In this section we 
prove the following theorem. 

Theorem 7.1. Assume thai /io G R\cr(Aj2) is not an eigenvalue of the selfadjoint 
Dirichlet Laplacian A in L'^{Q), the scaling parameter A G T>a is not real, and 

Q<I3 < TO«n^e<,(An)|3{(l + >^)^/JI^~^}\■ (7.1) 

Let s(^, rf) stand for the scaling function IR.+ x 9 (a;, y) i— >■ s{x) written in the coor- 
dinates {C,ri) G C and extended to G by zero; see (2.2)-(2.4). Then for all sufficiently 
large r > and all F G L'^{Q) satisfying eP^T G i^(^J) the estimate 

\\ef^\Ax,r-li^)-^T;Hl{Q)\\ < C\\e^^J^;L^g)\\ (7.2) 
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is valid with a constant C independent of T. 

Theorem 7.1 together with Theorem 6.1 shows that under the additional assump- 
tion e<^\! o -dx^r) e L'^{Q) on f e UaiG) infinite PMLs absorb out going or mcommg 
solutions (depending on the sign of 9A) so effectively that the function UA,r in Theo- 
rem 7.1.3 exponentially decays at infinity in the mean. 

Proof. Consider the conjugated operator e'^^A^.re"'^^ as an unbounded operator 
in L'^{G) with the domain C§°{Q). With this operator we associate the quadratic form 
qxrW'ju] = {e^^Ax^re~^^u,u)x^r- Observe that 

^[w,-u]-gA,r[w,w] = -/3^ / {exl-i^V(r,s,uV(^ns)dCdri 
Jg 

-0 ((deteA,r-)^^^e^^^uV^^s, Vf,,(detei^ J"^^\)dC(ir? 
Jg 

+ /5 / {^xl-'^Cv'^''^'^(ris) dC dr], 
Jg 

where qx,r is the same as in (4.1). Since the right hand side depends linearly on Vf,jU, 
similarly to the the proof of Lemma 4.3 one can deduce 

\(lx,MM - qx,r[uM\ < e|(-Vc^ • e^_i^V^^u,u)| + Ce-^u; L''{g)f 

with an arbitrary small e > and a constant C independent of w G C^{g). This 
estimate together with (4.8), (4.9), and Lemma 4.1 implies that for all sufficiently 
large r > we have 

\&rg{ql^[u,u]+'y\\u;L\g)f)\<^ (7.3) 

with some angle ip < n/2 and 7 > 0, which are independent of u e C^{Q). Therefore 
e^^Ax^rG^'^^ with the domain C§°{g) is a densely defined sectorial operator in L'^{Q). 
Let D(e'^^AA,re~^^) be the domain of its m-sectorial Friedrichs extension [23, Chapter 
VI.2]. The Friedrichs extension will also be denoted by e'^^AA.re"''^. 

As in the proof of Proposition 4.4.1 we conclude that C^{Q) is a core of the 
m-sectorial operator e^^AA,re~'^^. In Section 8 we will show that the graph norm of 
e'^'AA.r-e"'^' is equivalent to the norm in Hl{g); hence D(e'^=AA,re"'^') = Hl{g). Fur- 
thermore, we will localize the essential spcctnmi ae,ss{pJ^^ Ax.r^^ ^^^) of the m-sectorial 
operator e^^AA,re~^^, see Proposition 8.1. It turns out that the essential spectrum 
consists of an infinite number of parabolas, see Figure 7.1. In the case /3 = the 
parabolas collapse to the dashed rays originating from the thresholds v G cr(An) and 
we obtain the essential spectrum aess{Ax,r)- 

The m-sectorial operator e'^^AA.rC"''^ defines the Predholm holomorphic operator 
function 

M ^ e^'Ax,re-^' - M : H^{g) ^ L\g) 

on the simply connected subset of C \ aess{e^^Ax,re^^^) containing an infinite part 
of the real negative semiaxis (regular points of e^^^AA^re '^^)- Condition (7.1) on ^ 
guarantees that the point /io is in this simply connected subset. As the spectrum 
of a Fredholm holomorphic operator function consists of isolated eigenvalues of finite 
multiplicity, /ig is a regular point or an eigenvalue of e^^AA^rfs^'^^. The inclusion 
* e ker(e'^^AA,re"'^^ - /Uo) implies e"'^^* e ker(AA,r - A*o), and hence * = by 
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Fig. 7.1. Essential spectrum of the conjugated operator e^^Ax^j.e for OA > and /9 ^ 0. 



Theorem 6.1.1. Thus the operator e^^Ax.rS — f^o yields an isomorphism between 
the spaces Hq{Q) and L^{G)- This together with the equahty 

justifies the estimate (7.2). □ 

8. Localization of the essential spectrum. In this section we localize the 
essential spectrum aess{&^^^\,re~^^) of the m-sectorial operator e^^Ax,re~^^ in L'^iQ) 
with parameters A e and /3 e M; here s is the same as in Theorem 7.1. In 
particular, in the case = wc find (Tess{Ax,r). We also prove that Die^^Ax.r-e^'^^) = 
Hq{Q). In other words, we show that the information on the essential spectrum and 
the domain of e'^^A^^^e"''^ we used in the proofs of Theorem 6.1 and Theorem 7.1 is 
correct. 

Let us note that for fixed A and /3 the spectrum aess{e'^^^x,re~^^) depends only 
on the behavior of the scaling function s and the matrix ex,r outside any compact 
region of Q. In order to control (7ess{p-^^Ax.re~^^) we imposed the condition (2.4). 

Proposition 8.1. Assume that A e Va, (3 G M., and r > is sufficiently large. 
Then the following assertions hold. 

1. The Hilbert spaces D{e^^Ax,re~^^) and Hq{Q) are coincident and their norms 
are equivalent. 

2. The hounded operator 

e^'Ax,re-^' - /i : H^iG) ^ L^G) (8.1) 

is not Fredholm (or, equivalently, ji S aes8{p^^^\,re-~^^) ) if and only if the 
parameters fjb, \, and /3 meet the condition 

- // = (1 + A)"2(/3 + i£f for some v e cr(An) and C G M. (8.2) 

Proof. The proof is essentially based on methods of the theory of elliptic non- 
homogeneous boundary value problems [24, 25, 28, 26]. We will rely on the following 
lemma due to Peetre, see e.g. [26, Lemma 5.1], [25, Lemma 3.4.1] or [32]: 

Let X,y and Z he Banach spaces, where X is compactly emhedded into Z. 

Furthemore, let C he a linear hounded operator from X to y. Then the next 

two assertions are equivalent: (i) the range of C is closed iny anrfdimker^ < 

cxD, (ii) there exists a constant C , such that 



\\u;X\\<C{\\Cu;y\\ + \\u;Z\\) Vw e A". 



(8.3) 
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Below we assume that /x, A, and /3 does not meet the condition (8.2) and establish 
the coercive estimate 

\\u;H^ig)\\ < C{\\{e^'Ax,re-^'-^l)u■,L^{g)\\ + \\^u;L\g)\\) \/u e H^iG). (8.4) 

Here w e C°°(5) is a positive rapidly decreasing at infinity weight, such that the 
embedding of Hq{Q) into the weighted space L^{Q;w) with the norm ||w-; Z/^(^)|| is 
compact. Note that (8.4) is an estimate of type (8.3) for the operator (8.1). 

The strongly elliptic differential operator e^^Ax^r^^'^^ endowed with the Dirichlet 
boundary condition set up a regular elliptic boundary value problem. Solutions of 
a regular elUptic boundary value problem satisfy local coercive estimates, e.g. [26] 
or [25]. Thus we have the local coercive estimate 

Wptu; HUg)\\ < C{\\gT{e^' Ax,re-^' - fi)u; L^g)\\ + Wqtu; L^{g)\\). (8.5) 

Here px and qt arc smooth compactly supported cutoff functions in Q such that 
Pt(C)??) = 1 for 1(1 < T + 1 and QtPt = pr, where T is a large fixed number. 

Let XT e C°°{g) be another cutoff function such that XTiC^v) = 1 for \(\ > T 
and XriCiV) = for \(\ < T — 1. On the next step we establish the estimate (8.4) 
with u replaced by xtu. We will do it in the coordinates {x, y) G IR+ x ^. 

Let L^(R X il) be the space of functions in the infinite cylinder R x f2 with the 
norm (Jj^ ||u(x); Z/^(0)||^ dx) . Introduce the Sobolev space Hq{M. x f2) of functions 
with zero Dirichlet data on R x dQ, as the completion of the set Co°(R x O) with 
respect to the norm 

||u;ifo'(Kxfi)|| = ( \\did^u;L\^xQ)ff'\ 

t+\m\<2 

Denote u = {xtu) ° where h is the diffeomorphism (2.1). Let 

AA,r = -(detgA,r) V^jy • (dct gA,r) ^^^g^r^^y 7 ASPa, (8.6) 

be the operator A\^r written in the coordinates {x,y). Here g\^r is the matrix (3.3) 
and ~ {dx,dy^ . . -dy^)^ . Due to our assumptions on x the estimates < e < 
det V(x,7/) < 1/e hold uniformly in {x,y) e R+ x f2. Hence for some independent of 
u e Cq^{Q) constants ci, C2, and C3 we have 

\\XTu;Hl{g)\\ = \\A{xTu);L\g)\\ + \\xTu;L\g)\\ 

< ci(||Ao,rU;i^(R X f2)|| + ||u;L^(R x f2)||) < C2\\u;HI{^ x (8.7) 
\\{e^'Ax,re-^' - iJ,y,L\R x < C3||(e'3=AA,.e-''= - /x)xt«; L\g)\\. 

Here the functions u, s, and AA,rU = (AA,r(XTw)) o x are extended from x Cl 
to the infinite cylinder R x II by zero, and |lAo,rU;L^(R x < C||u;i?^(R x fi)|| 
because the coefficients of the Laplacian Ao,r are bounded, cf. (8.6) and (3.5). As T 
is large, the function u is supported in a small neighborhood of infinity. Due to the 
stabilization condition (3.4) on g^^ and the condition (2.4) on the scaling function s 
the coefficients of the differential operator 



e^'Ax,re-^' - An + (1 + X)-^{d, + 
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are small on the support of u. As a result we get the estimate 

II (e^=AA,.e-'3= - An + (1 + Xy^O^ + l3)^)u; L^R x < e||u; H^{R x (8.8) 

where e is small and independent of u e Cq^{G); moreover, e — >■ as T — )• +oo. 
Consider the bounded operator 

An - (1 + A)-2(a^ +Pf - n: H^{Rxn) ^ L^{R x n). (8.9) 

Applying the Fourier transform we pass from the operator (8.9) to the selfad- 

joint Dirichlet Laplacian Afj + (1 + A)~^(/3 + i^)^ — /x in L^(f2). Since /z, A, and /3 do 
not meet the condition (8.2), the spectral parameter /i — (1 + A)~^(/3 + z^)^ is outside 
of the spectrum of An for all ,^ e M. Then a known argument [25, Theorem 5.2.2], [24, 
Theorem 2.4.1], which is also used as a part of the proof of Lemma 9.2 below, implies 
that the operator (8.9) realizes an isomorphism. In particular, the estimate 

]]u; Bli^ X n)]] < c|| (An - (1 + A)-2(a^ + - m)u; ^^(M x fi) || 

is valid with an independent of u S ifo(^ ^ ^) constant c. As a consequence of this 
estimate and (8.8) we obtain 

(1 - ec)]]u; Jf2(R x O)]] < c|| (An - (1 + A)-2(a^ + - /x)u; ^^(M x 0)|| 
- c|| [e^' l\x,re-^' - An + (1 + A)-2(a^ + /3)2)u; ^^(M x 0)|| 
< clKe^^AA.r-e"'^" - m)u;L^(K x fi)]]. 

If T is sufficiently large, then ec < 1. This together with (8.7) gives 

\\XTU-,Hlm<C\\{e^'Ay^,rer^^ - ^i)XTU-,L\g)t (8.10) 

where the constant C = c(l — ec)~^C2C3 is independent of u e Cq°(5). By continuity 
the estimate (8.10) extends to all u € Hq{Q). 

Now we combine (8.10) with (8.5), and arrive at the estimates 

\\u;H^{g)\\ < \\xTu;H^{g)\\ + \\ptu; H^{g)\\ 
< C(llxT(e'5=AA,.e-'5=-M)«;i'(a)|| + II [e'^=AA,.e-'5=,XT]«; 1.2(^)11 
+\\QT{e^'A^,re-^' - /x)w; L\g)\\ + Wqtu; L\g)\\) 
< C{\\{e^'Ax,re-^'-ti)u;L\g)\\ + \\gTu;L\g)\\). 

Here we used that = 1 on the support of the commutator [e^^Ax^re~^^,XT], and 
hence 

\\[e^'Ax,re-^^XT]u;L\g)\\ < C\\ptu; H^{g)\\. 

For an arbitrary positive weight w wc have H^Tw; -^^(5)|| < C||wu; L^(^)|| with an 
independent of m S -ffp (^) constant C. Thus the estimate (8.3) is a direct consequence 
of (8.11). By the Peetre's lemma we conclude that the range of the operator (8.1) is 
closed and the kernel is finite-dimensional. 

Clearly, the graph norm \\u; L^{g)\\ + \\el^'Ax^re-'^'u\\ of u G C(^{g) is majorized 
by \\u;HQ{g)\\. The estimate (8.3) with w = 1 implies that the norm iJg (C/)|| is 
majorized by the graph norm of u. Since the set Co°(^) is dense in D(e'^^AA,re~''^) 
and in HQ{g), this proves assertion 1. 
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In order to see that the cokernel coker{e^^Ax^r^~^^ — i^) = ker[{e^^Ax^rS~^^)* —Jl) 
of the operator (8.1) is finitc-dimcnsional (if /x, A, and /3 does not meet the condi- 
tion (8.2)) we derive the coercive estimate 

\\u;H^{g)\\ < C{\\{{e^'Ax,re-^r-Ti)u;L^{Q)\\ + \\v^u; L\g)\\) (8.12) 

for the adjoint {e^^A\^re~^^)* of the m-sectorial operator e^^A\^re-~^^ and apply tlic 
Peetre's lemma. The m-sectorial operator (e'^^A^^re^'^^)* corresponds to the closed 

densely defined sectorial form qx^r[u,u] with the domain H^{Q). The proof of the 
estimate (8.12) is similar to the proof of (8.4), we omit it. 

Wc have proved that the operator (8.1) is Frcdholm provided the condition (8.2) 
is not satisfied. Now we assume that the condition (8.2) is met, and show that the 
operator (8.1) is not Fredholm. 

Let X be a smooth cutoff function on the real line, such that xi^) = 1 for l^^— 3| < 1 
and x{x) = for |a; — 3| > 2. Consider the functions 

ue{x,y)=x{x/e)exp{i{l + X)y/jr^x-l3x)^y), {x,y)eRxn, (8.13) 

where $ is an eigenfunction of An, corresponding to the eigenvalue v € cr(An). The 
exponent in (8.13) is an oscillating function of x. Straightforward calculation shows 
that 

\\{Aa-{l + X)-^{d^ + ^f-lj)ue;L''{Rxn)\\<C, ||u^; ir^(R x 1^)|| ^ oo (8.14) 

as ^ ^ +00. Similarly to (8.8) we conclude that 

\\{e^'Ax,re-^'-Aa + {l + X)-\d,+l3f)ue;L\Rxn)\\<€(\\ue;H^{Rxn)\\, (8.15) 

where ^ as £ — >■ +00. Let the functions U£ = U£ok~^ be extended from C to ^ by 
zero. If, on the contrary, the operator (8.1) is Fredholm, then by the Peetre's lemma 
the estimate (8.4) holds with any weight w, such that Hq{Q) ^ L'^{Q] w) is a compact 
embedding. Without loss of generality we can assume that ||wu£; L^(t/)|| < C for all 
i>l. After the change of variables {Cv) {^jV) the estimate (8.4) implies 

||u^; H^{R X n)\\ < C{\\{e^'Ax,re-'^' - i^W,L\R x f2)|| + 1), 

where the function e^''Ax,r^^'^'^^^ = (e'^^^A.re '^^u^) o >c is extended from M+ x 51 to 
M X f2 by zero. This together with (8.15) justifies the estimate 

II u,; X fi)|| < C{\\{An - (1 + X)-^{d^ + - A<)u^;i'(K x r!)|| + l), 

which contradicts (8.14). □ 

9. Problem with finite PMLs. Consider the truncated domain Qr with piece- 
wise smooth boundary, see (2.9). Introduce the Sobolev space H^{Qii) as the comple- 
tion of the set C^iGR) in the norm \\v:Hl(gR)\\ - {Y.i+\m\<2 ; L2(g^)||2)i/2, 
In this section wc study the problem with finite PMLs: Given g € L^i^Qn) find a 
solution V e Hq{Qr} of the equation 

{Ax,r - l^o)v = g in Qr, (9.1) 

where R> r. The next theorem presents a stability result for this problem. 
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Theorem 9.1. Assume that Ho S R\(j(Aq) is not an eigenvalue of the selfadjoint 
Dirichlet Laplacian A in L'^{Q). Take a sufficiently large r > and A G \R. Then 
there exists Rq > r such that for all R> Rq and g G L^{Gr) the equation (9.1) has a 
unique solution v € Hq{Qb)- Moreover, the estimate 

\\vMl{gn)\\<C\\g-L\gR)\\ (9.2) 

holds with an independent of R> Rq and g constant C. 

The proof of Theorem 9.1 will be carried out by the compound expansion tech- 
nique. This requires construction of an approximate solution to the equation (9.1) 
compounded of solutions to limit problems. As the first limit problem we take the 
problem with infinite PMLs. As the second limit problem we take a Dirichlet problem 
in the scmi-cylindcr (— oo, R) x CI studied in the next lemma. 

Lemma 9.2. Introduce the weighted Sobolev space ifg^^ ((— oo, R) x f2) of functions 
satisfying the homogeneous Dirichlet boundary condition as the completion of the set 
Cq°{{—oo,R] X Q) with respect to the norm 

\\u;Hl^{{-oo,R) X 0)11 = ( ^ T \\e-P^did^u{xy, L'{n)f dx 

\e+\m\<2-'-'^ 

Let I/^((— oo, i?) X O) be the weighted L"^ -space with the norm 

||f;L|((-oo,i?) X 0)11 = ( /"^ ||e-'^^f(a;);L2(0)||2dx 

Assume that X G Va \ M, /kq G K \ it(Ao), and j3 is in the interval (7.1). Then for 
any f G L|((— oo,i?) x O) there exists a unique solution u G Hqq(^{—oo,R) x O) to 
the equation 

(Aa-(l + A)-2a^-Mo)u = f. (9.3) 

Moreover, the estimate 

\\u;Hli3{i~oo,R) X 0)11 < C\\f;Lj{{-oo,R) x 0)|| (9.4) 

holds, where the constant C is independent off and R. 

Proof It suffices to prove the assertion for R = 0. Then the general case can be 
obtained by the change of variables x x — R. 

The set C^((— oo,0) x O) of smooth functions with compact supports is dense 
in L^((— oo.O) X O). We first assume that f G C^((— oo,0) x O) and extend f 
to a function in C^(]R x O) by setting f(— x) = — f(.x) for x < 0. Consider the 
equation (9.3) in the infinite cylinder M x O. As is known [36], the Fourier transform 
f(^) = e"^f (a;) da; is an entire function of ^ with values in L^(0); it is rapidly 
decaying at infinity in any strip G C : < /?} in the sense that the estimates 
||f(^);L^(0)|| < C/3,ft(H-|^|)~'' hold for /c = 0, 1, 2, . . . and some constants C;3,fe. Since 
/3 is in the interval (7.1), the distance d between the set {/xo — (1+A)~^^^ : < < /3} 
and the spectrum ^(An) of the selfadjoint operator Aq in L^(0) with domain Hq{Q,) 
is positive. Hence for 

*(0 = (An + (1 + A)-^^' - Mo) "'f (0. < 3^ < /3, 





DIRICHLET LAPLACIANS IN QUASI-CYLINDRICAL DOMAINS 



23 



we have the estimate ||^'(^); L2(rj)||2 < d-'^\\f{^);L'^{n)\\'^. This together with the 
eUiptic coercive estimate 



for the Dirichlct Laplacian Aq gives 

mo-,H^m'<ic+d-')\m-,L'm^ o<3^<^. (9.5) 

The differential operator An — (1 + X)~^d^ is strongly elliptic. Therefore the local 
coercive estimate 

\\eu;H^{R X n)\\^ < c(||?f;X2(M x 17)f + ||?u;L2(m x (9.6) 

is valid, where g and <; arc smooth functions of the variable x with compact supports 
such that g ^ and gq = g. Wc substitute u{x,y) = e*^^\E'(^, t/) into (9.6). After 
simple manipulations we arrive at the estimate 

E l?l''l|5r*(0;i'(^^)lP<c(||f(0;i'(^^)f + Il*(0;i'(f2)ll'), (9-7) 

e+\m\<2 

where the constant C depends on g and but not on ^ or f . If |^| > C with sufficiently 
large C > 0, then the last term in (9.7) can be neglected. This together with (9.5) 
justifies the estimate 

e+\m\<2 

where < 3^ < /? and the constant C is independent of ^ and 'I'(^). Therefore the 
analytic in strip < 3^ < 3 function ^ h- > v]/(^) ^ (fi) is rapidly decaying at infinity. 
This together with the Cauchy integral theorem allows us to replace the contour of 
integration in the inverse Fourier transformation u(a;) = (27r)~^ J^e~'^^^^{^) d^. We 
obtain 

u{x) ={2n)-' [ e-'«-(An + (l + A)-2^2-Mo)"'f(Orfe 

The Parseval equality gives 

2n [ \\e-f'^dX<^y,L\m'dx= [ l^'lK^iO; L\m' d^, 

27r / \\e-^^f{x);L\n)fdx= [ ||f (0; L'(fi)f rf^- 
Integrating (9.8) with respect to ^, ^ — i/3 € M, we deduce the estimate 

/oo POO 
We-^'^did^uix); L'^{n)f dx<C ||e-''^f(a;); L'^{n)f dx; (9.9) 
. .. -OO J —OO 
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in the case /3 = this estimate takes the form II u;i/^(Mxfl) II < C\\f:L^{Rxn)\\. Thus 
for any f G C^(R x fl) there exists a solution u € Hq^R x fl) to the equation (9.3) 
and the estimate (9.9) holds with any /3 in the interval (7.1). Usual argument on 
smoothness of solutions to elliptic problems gives u G C°°{R x fl). Prom the equality 
f(-x) -f(x) it follows that f(-0 = -f(0 and therefore *(0 = -*(-C)- Hence 
u(a;) = —u{—x) and u(0) = 0. As in the proof of Proposition 4.4 we conclude that in 
the norm of Hqj^(^{0, oo) x fl) one can approximate u by functions in C^((0, oo) x il). 
Hence u e i?Q^((— oo,0) x il). By continuity our construction extends to all f G 
L?j((— oo.O) X rj). In particular, for any / G Lq((— oo,0) x £7) we can find a solution 
u e i?o,o((~c», 0) X fl) to the equation (9.3). If f € L^((— oo,0) x f2) with some 
/3 in the interval (7.1), then the estimate (9.4) with = is a direct consequence 
of (9.9). It remains to note that u G Hq q(^{—oo,Q) x Jl) is a unique solution as 
our argument also shows that for any f G Lq((— oo,0) x fi) the adjoint equation 
(Ao — (1 + A)~^9^ — /io)u = f is solvable in the space i?Q o((— oo, 0) x il). □ 
Now we are in position to prove Theorem 9.1. 

Proof. The scheme of the proof is similar to the one we used in the proof of [20, 
Theorem 4.1]. We rely on a modification of the compound expansion method [27]. 
We say that w = w{R) G H^{Gr) is an approximate solution of the problem with 
finite PML if the following conditions are satisfied: 

i. The estimate \\w] HQ{Qji)\\ < c\\g; L'^{Qfi,)\\ holds with an independent of g 
and R constant c; 

ii. The estimate \\{Ax,r — ^^o)w — g; L'^{Gr)\\ < CB]\g'- L^{Qr)\\ is valid, where 
the constant C'n is independent of g and Cn — >■ as i? — >■ +oo. 

Due to condition i w continuously depends on g. Condition ii implies that the discrep- 
ancy, left by w in the equation (9.1), tends to zero as i? — > +00. Once an approximate 
solution w is found, it is not hard to verify the assertion of the theorem. 

Let p G C°°(M) be a cutoff function such that p{x) = 1 for x < and xi^) ~ 
for X > 1/2. We set pn = p{x — R), gR{C,v) = PR ° x~^{C,v) iC^v) € C, and 
Qr{C,v) = 1 for iCv) G\C. Let J" = gR/2g and f = (5 - /) Igr We 
extend from Gr to Q and f from (0,i?) x to (— oo,i?) x by zero. We find 
an approximate solution w compounded of UA.r = {'^x,r — Mo)~^-^ and a solution 
u € ilg q((— 00, R) X n) to the equation (9.3) in the form 

W = QRUx.r + (1 - to/3)(u o 

here the second term in the right hand side is extended from Qr fl C to Qr by zero. 

Let us show that w is an approximate solution. Observe that on the support of 
/ we have e^^ < Ce^^/^ and on the support of f we have e~^^ < Ce~^^/'^ uniformly 
in R. Hence 

||e''^.F;L2(g)||+e^«||f;L2((_^,i?) xO)|| <Ce^«/2||g;L2(g)|| (9.10) 
with an independent of R and g constant C. Similarly to (8.7) we conclude that 
\\w; Hl{gR)f < WQRUxy, Hl{g)f + c(|| A;,,,(l - pR,^)yi- L^{{-oo, R) x O) || 
+ ||u;L2((_oo,i?) X n)\\y < C\\uxy,H^ig)f + C\\u;Hlo{(~^,0) x 0)f , 

where C is independent of R and A^.r is the operator (8.6). This together with the 
estimates (9.4), (9.10) for ^ = and Theorem 6.1.1 implies that the condition i is 
satisfied. 
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Let us verify the condition ii. We have 

(AA,r - Mo)w - 3 = [Ax^r, gR]u\,r + (1 + A)"^ ([9^, Pi?/3]u) o 



-((A,^, - A^, + (1 + Xr'dl) (1 - p«/3)u) o 



(9.11) 



The support of the term [A^.r, efl]MA,r is a subset of the image of {R,R + 1/2) x f2 
under the diffeomorphism x. On this support the weight e^^ is bounded from below 
by ce^^ uniformly in R> 0. As a consequence we get the uniform in R estimates 

\\[Ax,r, eR]uxy,L''{gn)\\ < C^e-^V^xy,H^iG)\\ < C2e-0Vj^;L^ig)\\, 

(9.12) 

where we used Theorem 7.1. Now we estimate the second term in the right hand side 
of (9.11). On the support of [dl,pii/3]u we have e"''^ > Ce'^^/^. Relying on (9.4) 

we obtain 

II (1 + X)-^{[dl pi,,/,]u) o ^-^:L\gn)\\ < c\\ [dl p^/glu; i'((-oo, R) x n) \\ 
<CieP''/^u;Hl^{{-oo,R)xn)\\<C2e^''/^f;L^{{-oo,R)xn)\\. ^' ^ 

FinaUy, consider the last term in the right hand side of (9.11). On the support of 
(1 — Pii/3)u the coefficients of the operator Ax,r — Aq + (1 + X)~^d^ tend to zero as 
R — )• +00, cf. (3.4) and (8.6). This together with the estimate (9.4) for ^ = gives 



((Aa,, - An + (1 + Xr^dl){l - p«/3)u) o >,-';L\gn] 
< CR\\f;Ll{{-oo,R) X fl)\\, 



(9.14) 



where ^ as i? — )■ +00. Prom (9.10)-(9.14) it follows that w meets the condition 
ii. Thus w = w{R) is indeed an approximate solution to the problem with finite PML. 

Now we are in position to prove the assertion of the theorem. Observe that 
{A\^r — ^io)w — g = 0{R)g with some operator D{R) in L'^{Qi{), whose norm |||D(i?)||| 
tends to zero as i? +00 because of the condition ii on w. For all i? > i?o with 
a sufficiently large i?o we have |||D(-R)||| < |||D(i?o)||| < 1- Hence there exists the 
inverse (/ + D{R))~^ : L'^(Qr) L'^{Gb) and its norm is bounded by the constant 
1/(1 - |||D(i?o)|||) imiformly in R > Rq. We set g = (I + 0{R)y^g. In the same way 
as before we construct the approximate solution w for the problem (9.1), where g is 
replaced by g. Then for v = w we have (A;^^^ — Mo)^^ = 9 + ^(^)5 = 9 and 

\\v;H^igR)\\ < crg;L\gR)\\ < c/{l - mR^)\\\)\\g; L^g^)]], 

where C is independent of > Rq. Thus for R > Rq and g G L'^{gii) there exists 
a solution v G H^{gfj) to the equation (9.1) satisfying the estimate (9.2), where the 
constant C is independent of R and g. In the remaining part of the proof we show 
that this solution is unique. 

Let A;^^ be the unbounded operator in L^(5i?) such that for any v in its domain 
HQ^gR) we have A^^v = Ax^rV. Note that A^^ is the operator of the problem with 
finite PML. To the operator A^ ^ there corresponds the quadratic form 

Qx,r['"''"^ = / {i<ietex,ry^^ex].VcnV,Vcn{dete:^^)~^^'^v)dCdT] 

in L'^{gf{) with the domain HQ{gfi). In the same way as in Section 4 we con- 
clude that the form q^^ admits a densely defined sectorial closure with the do- 
main H^{gii), where H^{gji) is the completion of HQ^gn) with respect to the norm 
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\\v;H\gR)\\ = {Ei+\m\<i\\9cd:;'v;L^{GRW)'/^. This together with the argument 
above impHcs that ah ^ < with sufficiently largo absolute value are regular points 
of the operator A^^. Hence the sectorial operator A^^ coincides with its m-sectorial 
Priedrichs extension. Moreover, thanks to the argument above we know that under 
the assumptions of theorem for any g G L'^{Gr) there exists v G Hq{Qr) such that 
(A;^^ — iX{))v = g. Similarly, one can study the adjoint m-sectorial operator (Aj^^)*. 
It turns out that under the assumptions of theorem for any g g L'^^Qr) there exists 

V in the domain H^iGR) of {Af jj* such that ((Af^)* - ij,o)v = g. Therefore the 
equation (9.1) is uniquely solvable in Hq{Qr). □ 

In the next theorem we show that under some natural assumptions solutions 

V = v{R) of the problem with finite PMLs converge in the domain Qr to outgoing 
or incoming solutions u± with an exponential rate as i? — > +oo. In other words, we 
estimate the error produced by truncation of infinite PMLs. Solutions to the problem 
with finite PMLs can be found numerically with the help of finite element solvers; 
certainly, discretization produces yet another error that we do not estimate here. 

Theorem 9.3. Assume that fiQ S R\a(An) is not an eigenvalue of the selfadjoint 
Dirichlet Laplacian A in L'^{G), the parameter r > is sufficiently large in the sense 
of Remark Jf..2, A G X'q \ R, and j3 is in the interval (7.1). Let f € 'Ha{G) satisfy 
the inclusion e^^{f o i?A,r) S L^{Q), where s is the same function as in Theorem 7.1. 
In (9.1) we set g = f o i?A,r- Then there exists Rq > r such that for R> Ro a unique 
solution V = v{R) G Hq(Qr) of the problem with finite PMLs converges in Qr 

1. to the outgoing solution U- € -f^o loc(^) ^/ equation (A — no)u = f in the 
case QX > 

2. to the incoming solution u+ e Hq i^^{Q) of the equation (A — /io)w = f in the 
case SjA < 

in the sense that as R^ +oo the estimate 

^ \\dld^{u± - vn); L^gr)r < Ce-^^^'We^^f o dx,r); L'ig)f (9.15) 

e+\m\<2 

holds with a constant C independent of R> Ro and f. 

Let us remark here that the assumptions of Theorem 9.3 on the right hand side / 
are a priori met for all / S L'^[Q) such that / \ C = F o k with some F S ; here <S is 
the algebra defined in Section 5 and >c is the diffeomorphism (2.1). Prom Lemma 5.2 
it follows that the set of functions / satisfying the assumptions of Theorem 9.3 is 
dense in L^{Q). In particular, we can take any / e L^{Q) supported in Qr. 

Proof. Let p G C°°{R) be a cutoff function such that p{x) = 1 for a; < and 
Xix) = for X > 1/2. We set pR = p{x - R), gniC, rj) = pro K~'^{C,rj) for (C, rj) e C, 
and gR{C,ri) = 1 for {C,ri) G G\C. Thanks to Theorem 6.1.3 it suffices to prove the 
estimate (9.15) with u± replaced by QRUx^r- The difference QRUx,r — vr G Hq{Qr) 
satisfies the problem (9.1) with g = {qr — 1)(/ o t?A,r) + [AA,r, ^ij]wA,r- Observe that 

\\{QR-l){fo^x,r);L\gR)\\<Ce-^''\\e^'{fo^^^r)-L\g)\\, 
mx,r,QR]uxy,L\gR)\\ < Ce-^''\\e^^uxy,Hl{g)\l 

because the functions [qr — 1)(/ o '&x,r) and [A^^r, QR\Ux,r^ being written in the co- 
ordinates {x,y), are equal to zero for X < R, while e^^^'i^ = ce^^^ for a; > ii > C, 
cf. (2.4). This together with Theorem 7.1 gives 



Wg; L\gR)\\ < ce-'3«||e'3=(/ o ^,,,); L\g)\\. 



(9.16) 
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By Theorem 9.1 we have 

\\0RUx,r-vii;H^{gn)\\<C\\g;L^{gn)\\, R>Ro. (9.17) 
It remains to note that 

^ ||a^a™(n± - VRy,L\gr)f < \\QRUx,r - VR;H^{gR)f, R>Ro>r. 

e+\m\<2 

This together with (9.17) and (9.16) completes the proof of the estimate (9.15). □ 

Remark 9.4. By Theorem 9.3 the rate of convergence of the PML method de- 
pends only on the spectral parameter /xq and the infinitely distant cross-section Q.. In 
the case C = (0,oo) x Q, the quasi- cylindrical domain g corresponds to a resonator 
with attached tubular waveguide. In this particular case the results of Theorem 9.3 
for f e L'^{g) with supp/ C t/r can equivalently be obtained by the modal expansions 
technique, e.g. [5, 6, 7]. 
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